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1  introduction 


Background 

A  large  number  of  U.S.  Army  training  and  testing  installations  are  located  in 
arid  and  semi-arid  environments  in  the  desert  Southwestern  United  States. 
Fort  Bliss,  Texas,  is  a  Training  and  Doctrine  Command  (TRADOC)  installation 
located  in  the  Northern  Chihuahuan  Desert  of  Western  Texas  and  South  Central 
New  Mexico.  At  1.1  million  acres,  it  represents  the  single  largest  TRADOC  in¬ 
stallation.  Fort  Bliss  exists  within  a  semi-arid  ecosystem;  therefore  it  is  highly 
susceptible  to  disturbance  caused  by  both  natural  and  anthropogenic  factors. 
Disturbance  in  arid  landscapes  is  usually  associated  with  a  decrease  in  vegeta¬ 
tive  cover  and  an  increase  in  soil  erosion,  and  can  occur  at  many  different  scales. 
Military  training  (both  tracked  and  wheeled),  grazing,  and  recreational  activities 
all  act  as  stressors  on  the  landscape  at  Fort  Bliss.  The  Fort  Bliss  Directorate  of 
Environment  (DOE)  is  responsible  for  sustainment  of  natural  resources  in  sup¬ 
port  of  the  training  mission.  Therefore,  Fort  Bliss  natural  resource  managers 
require  information  on  land  condition  to  make  informed  land  management  deci¬ 
sions  and  to  support  conservation  and  compliance  efforts,  including  direct  sup¬ 
port  of  Environmental  Impact  Statements  (EISs)  for  withdrawal  of  Bureau  of 
Land  Management  (BLM)  lands  for  training.  Fort  Bliss  natural  resource  man¬ 
agers  require  a  cost-effective  method  for  assessing  and  monitoring  land  condi¬ 
tion. 

Vegetation  amount,  condition,  and  composition  are  the  most  important  indicators 
of  land  condition  in  arid  environments;  therefore,  it  is  necessary  to  conduct  in¬ 
ventories  and  monitor  vegetation  amount  and  condition  (Robinove  et  al.  1981; 
Musick  1984;  Duncan  et  al.  1993).  Many  Department  of  Defense  (DOD)  installa¬ 
tions  have  implemented  the  Land  Condition  Trend  Analysis  (LCTA)  program, 
which  is  part  of  the  Integrated  Training  Area  Management  (ITAM)  plan  for  the 
U.S.  Army.  LCTA  provides  a  standard  method  of  inventory  and  monitoring  for 
vegetation  and  wildlife  (Tazik  et  al.  1992).  Permanent  plots  are  established  and 
visited  annually  to  conduct  a  detailed  census  of  vegetation  and  wildlife.  How¬ 
ever,  field  surveys  are  expensive;  therefore,  a  complete  survey  of  large  installa¬ 
tions  is  cost  prohibitive  and  a  detailed  vegetation  map  is  often  lacking.  Long¬ 
term  trends  in  vegetation  condition  can  be  monitored  by  LCTA  field  surveys,  but 
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it  is  not  feasible  to  assess  vegetation  condition  at  any  one  time  over  a  large  area 
like  Fort  Bliss  based  solely  on  field  surveys. 

Remotely  sensed  imagery,  because  of  its  large  areal  coverage  and  high  temporal 
frequency,  provides  an  ideal  supplement  to  field  surveys  when  attempting  to 
characterize  and  monitor  changes  in  vegetation  condition  at  different  scales  and 
levels  of  detail.  Archival  coarse  resolution  imagery  can  be  used  to  detect  relative 
changes  in  surface  condition  or  vegetation  cover  in  both  small  and  large  areas. 
Vegetation  indices  and  soil  brightness  indices  calculated  from  remotely  sensed 
reflectance  data  are  sensitive  to  the  amount  of  vegetation  on  the  Earth’s  surface. 
Special  indices  have  been  customized  to  address  unique  problems  associated 
with  reflectance  characteristics  of  arid  environments  where  vegetation  cover 
typically  is  low.  Assessing  trends  in  vegetation  and  brightness  indices  calculated 
from  temporal  image  data  sets  provides  a  cost-effective  method  for  monitoring 
relative  trends  in  land  condition.  Quantification  of  land  condition  with  remotely 
sensed  imagery  is  difficult  and  requires  an  empirical  relationship  between  indi¬ 
ces  and  ground  measurements  for  the  purpose  of  calibration.  Historical  ground 
surveys  are  not  always  available.  However,  documentation  of  historical  trends 
in  land  condition  is  useful  to  upper  management  for  evaluating  the  effectiveness 
of  historical  resource  meinagement  efforts,  including  Land  Rehabilitation  and 
Maintenance  (LRAM)  over  long  periods  of  time  and  across  large  geographic  ar¬ 
eas.  Documentation  is  also  useful  for  generating  support  for  continued  efforts  or 
identifying  other  “problem”  areas  where  future  resource  management  efforts 
should  be  focused.  Detecting  coarse  resolution  changes  may  provide  general 
measures  of  significant  change  in  percent  cover  for  a  training  area  or  group  of 
training  areas,  for  example,  but  may  not  be  able  to  detect  more  subtle  changes  in 
grass  vs.  shrub  cover. 

Objective 

The  objective  of  this  research  was  to  characterize  the  small  scale,  gross  level 
change  in  land  condition  over  23  years  for  a  selected  area  of  Fort  Bliss,  Texas, 
based  upon  evaluation  and  analysis  of  several  vegetation  and  brightness  indices 
calculated  from  temporal,  archival  multispectral  imagery. 


Approach 

A  study  area  was  defined  to  include  most  of  Maneuver  Areas  5C,  5D,  and  parts  of 
4B  and  4C  on  the  west  side  of  U.S.  Route  54  and  an  adjoining  area  extending  ap¬ 
proximately  4  to  5  kilometers  (km)  east  across  Route  54  onto  the  McGregor 
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Guided  Missile  Range.  The  study  area  represents  a  gradient  of  impact  and  was 
selected  in  consultation  with  the  Fort  Bliss  Directorate  of  Environment  (Mr. 
Kevin  von  Finger).  The  maneuver  areas  have  been  impacted  to  varying  degrees 
during  the  time  period  of  this  study,  while  the  area  extending  onto  the  McGregor 
Range  has  been  relatively  unimpacted  by  tracked  and  wheeled  maneuvers  dur¬ 
ing  this  same  time  period.  An  archival  imagery  set  was  acquired,  georeferenced, 
and  subsetted  to  the  study  site.  Several  vegetation  and  brightness  indices  were 
reviewed  for  their  applicability  in  monitoring  land  condition  in  arid  environ¬ 
ments.  Based  on  a  literature  survey,  Red  Reflectance,  Modified  Soil-Adjusted 
Vegetation  Index  (MSAVI2),  Albedo,  and  the  Tasseled  Cap  Brightness  Index  were 
selected  for  analysis  and  calculated  for  19  image  dates  over  a  23-year  period. 
Trends  in  these  indices  were  evaluated  for  heavily  impacted  maneuver  areas,  low 
impact  areas,  and  the  entire  study  site.  These  trends  were  used  as  surrogate 
measures  of  vegetation  cover  and  land  condition  for  this  23-year  period  to  char¬ 
acterize  relative  changes  that  may  have  occurred. 


Scope 

The  specific  indices  and  resulting  trends  analyzed  in  this  report  are  specific  to 
Fort  Bliss,  Texas.  However,  assessment  of  relative  trends  in  land  condition 
based  on  vegetation  and  brightness  indices  calculated  from  temporal,  multispec- 
tral  imagery  is  a  valid  and  cost-effective  technique  that  should  be  applicable  to 
many,  but  not  all  military  installations.  These  techniques  work  best  in  environ¬ 
ments  with  a  relatively  open  canopy.  In  forested  ecosystems,  vegetation  indices 
may  provide  some  measure  of  land  condition  based  on  reflectance  of  the  upper 
most  strata  of  vegetation  or  forest  canopy,  but  would  not  be  sensitive  to  under¬ 
story  vegetation  or  vegetative  ground  cover  if  the  canopy  cover  exceeds  a  certain 
density.  The  temporal  period  of  23  years  (1972  through  1994)  in  this  study  rep¬ 
resents  nearly  the  maximum  amount  of  time  for  which  civilian  satellite  imagery 
has  been  available.  Archival  imagery  of  the  same  spatial  and  spectral  resolution 
should  be  available  for  most  DoD  training  and  testing  installations.  Therefore, 
the  techniques  applied  in  this  study  should  be  applicable  to  other  installations 
for  assessing  relative  trends  in  vegetation  cover  and  land  condition,  even  if  his¬ 
toric  ground  surveys  are  not  available.  The  same  techniques  can  be  applied  to 
future  imagery  to  continue  monitoring  land  condition. 

To  quantify  changes  in  land  condition,  a  correlation  between  spectral-based  indi¬ 
ces  and  ground-based  surveys  must  be  established.  In  a  recent  study,  correla¬ 
tions  were  established  between  several  vegetation  indices  and  various  LCTA 
groimd  surveys  (Zhuang,  Shapiro,  and  Bagley  1993;  Wu  and  Westervelt  1994; 
Senseman,  Bagley,  and  Tweddale  1996).  However,  due  to  the  lack  of  field  data 
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corresponding  to  the  imagery  used  in  this  project,  it  was  only  possible  to  assess 
relative  changes  in  vegetation  cover  and  land  condition  over  time. 


Mode  of  Technology  Transfer 

This  technical  report,  data,  results,  and  algorithms  used  represent  the  primary 
means  of  technology  transfer  to  Fort  Bliss  staff. 
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2  Background  on  Satellite  Data 

Common  commercial  satellite  sensors  record  reflectance  from  the  Earth’s  surface 
in  several  regions  of  the  electromagnetic  spectrum.  Depending  on  the  intended 
application,  different  wavelengths  are  best  suited  for  analyzing  different  aspects 
of  the  Earth-atmosphere  system.  For  the  purpose  of  assessing  vegetation,  the 
red  and  near  infrared  regions  of  the  spectrum  are  most  commonly  analyzed. 

Within  each  portion  of  the  spectrum,  different  properties  of  vegetation  control 
the  amount  of  electromagnetic  energy  that  is  absorbed,  transmitted,  or  reflected. 
In  general,  healthy  vegetation  has  been  characterized  by  low  reflectance  in  the 
visible  wavelengths  (400  to  700  micrometer  [pm])  and  high  reflectance  in  the 
near  infrared  wavelengths  (Kauth  et  al.  1978;  Tucker  1979;  Curran  1980).  Low 
reflectance  in  the  visible  wavelengths  is  due  primarily  to  high  absorption  of  pig¬ 
ments  such  as  chlorophyll  in  the  leaves.  Low  reflectance  in  the  blue  and  red 
bands  corresponds  to  the  two  chlorophyll  absorption  bands.  However,  there  is  a 
slight  increase  in  reflectivity  in  the  visible  green  portion  of  the  spectrum  (-550 
pm),  which  explains  why  healthy  vegetation  typically  appears  green  to  the  hu¬ 
man  eye. 

In  the  near  infrared  region  of  the  spectrum  (700  to  1300  pm),  cell  structure 
dominates  the  spectral  response.  Unlike  the  visible  regions,  there  is  no  strong 
absorption  in  these  wavelengths  and  they  show  relatively  high  reflectance  and 
scattering.  Therefore,  high  reflectance  in  the  near  infrared  wavelengths  is  di¬ 
rectly  proportional  to  plant  biomass.  In  general,  an  inverse  relationship  exists 
between  reflectance  in  the  visible  region,  particularly  in  the  red  wavelengths, 
and  biomass  production  of  a  plant,  while  reflectance  in  the  near  infrared  region 
of  the  spectrum  is  directly  proportional  to  plant  biomass  and  condition  (Campbell 
1987). 

Vegetation  Indices 

Several  vegetation  indexes  have  been  developed  to  reduce  multispectral  scanner 
data  observed  by  satellites  to  a  single  number  or  index,  for  the  purpose  of  quali¬ 
tatively  and  quantitatively  assessing  vegetation  conditions  (Tucker  1979;  Curran 
1980).  Almost  all  vegetation  indices  are  transformations  based  on  the  near  in¬ 
frared  and  red  portions  of  the  electromagnetic  spectrum.  In  general,  there  are 
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two  types  of  vegetation  indices:  ratio  and  orthogonal.  Ratio  indices  exploit  the 
contrasting  low  red  reflectance  and  high  near  infrared  reflectance  of  vegetation 
by  simple  ratios  of  these  two  bands.  Orthogonal  indices  are  based  upon  the 
“Tasseled  Cap”  transformation,  which  is  a  characteristic  plot  of  red  reflectance 
(x-axis)  vs.  near  infrared  reflectance  (y-axis)  that  is  useful  for  extracting  relative 
greenness  of  vegetation  and  soil  brightness  (Kauth  et  al.  1978).  Within  the  Tas¬ 
seled  Cap  distribution,  Kauth  determined  that  the  distribution  of  soil  reflectance 
variation  was  confined  to  a  “line  of  soils”  extending  from  the  origin  at  approxi¬ 
mately  45  degrees  from  the  x-axis  (red)  with  soil  brightness  increasing  with  dis¬ 
tance  from  the  origin.  Reflectance  variation  of  vegetation  is  then  measured  per¬ 
pendicularly  from  this  line  of  soils  in  the  direction  of  the  y-axis  (near  infrared). 
This  distribution  was  aptly  named  because  if  viewed  in  3  dimensions,  it  resem¬ 
bles  a  cap  with  a  tassel  extending  from  the  top.  Although  the  original  Tasseled 
Cap  transformation  was  based  upon  Landsat  Multispectral  Scanner  (MSS)  data, 
the  same  transformations  have  been  applied  to  Landsat  Thematic  Mapper  (TM) 
data  and  other  data  sets.  Regardless  of  the  sensor,  greenness  and  soil  brightness 
indices  derived  from  this  characteristic  plot  are  commonly  referred  to  as  Kauth- 
Thomas  or  Tasseled  Cap  Soil  Brightness  and  Greenness  indices. 

There  are  numerous  derivations  of  ratio  and  orthogonal  vegetation  indices  that 
are  beyond  the  scope  of  this  report.  However,  all  indices  are  similar  in  that  they 
provide  dimensionless  values  that  represent  relative  ranges  of  vegetation 
amount  or  condition.  In  general,  vegetation  indices  have  been  correlated  with  a 
number  of  vegetative  characteristics  such  as  biomass  (Tucker  1979),  percent 
cover  (Boyd  1986;  Senseman,  Bagley,  and  Tweddale  1996),  and  leaf  area  index 
(Richardson  and  Wiegand  1977). 

The  Normalized  Difference  Vegetation  Index  (NDVI)  is  probably  the  most  com¬ 
monly  applied  vegetation  index  for  assessing  vegetative  amount  and  condition 
(Rouse  et  al.  1973).  In  the  first  comprehensive  study  of  NDVI  and  vegetative  pa¬ 
rameters,  high  coefficients  of  determination  for  a  simple  linear  regression  were 
found  between  NDVI  and  total  wet  biomass,  total  dry  biomass,  leaf  water  con¬ 
tent,  dry  green  biomass,  and  total  chlorophyll  for  clipped  blue  grama  prairie 
grass  plots  (Tucker  1979).  NDVI  and  several  other  vegetation  indices  have  also 
been  correlated  with  abundance  of  arid  shrub  cover  in  northeast  Colorado  (An¬ 
derson,  Hansen,  and  Haas  1993),  Central  Texas  (McDaniel  and  Haas  1982),  Cen¬ 
tral  Washington  (Senseman,  Bagley,  and  Tweddale  1996)  and  specifically  in  the 
Chihuahuan  Desert,  where  NDVI  was  correlated  with  abundance  of  mesquite, 
creosote,  and  tarbush  (Duncan  et  al.  1993;  Musick  1984, 1986). 

However,  a  limitation  of  the  application  of  NDVI  in  arid  environments  is  that 
spectral  response  of  the  exposed  soils  often  dominates  the  spectral  response  of 
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vegetation.  This  is  due  to  sparse  vegetation  cover  and  high  near  infrared  reflec¬ 
tance  of  arid  soils  (Graetz  and  Gentle  1982).  Therefore,  several  indices  have 
been  developed  to  attempt  to  correct  for  this  factor,  including  the  Weighted  Dif¬ 
ference  Vegetation  Index  (WDVI,  Richardson  and  Wiegand  1977),  the  Soil- 
Adjusted  Vegetation  Index  (SAVI,  Huete  1988),  and  the  Modified  Soil-Adjusted 
Vegetation  Index  (MSAVI,  Qi  et  al.  1994).  These  indices  are  designed  to  maxi¬ 
mize  the  influence  of  vegetation  and  minimize  the  effect  of  background  soil. 
These  indices  require  a  user-defined  soil  correction  factor  “L”  as  input  into  the 
indices.  The  constant  L  represents  an  estimate  of  percent  vegetative  cover,  and 
therefore  is  often  unknown.  Typically,  different  L  factors  are  tested  and  the  final 
value  is  selected  based  on  agreement  with  ground  estimates,  or  an  L  factor  is  as¬ 
signed  based  on  user  knowledge  of  the  study  site.  In  general,  the  L  factor  is  dif¬ 
ficult  to  objectively  quantify.  The  Modified  Soil-Adjusted  Vegetation  Index  cal¬ 
culates  a  self-adjustable  L  factor  calculated  directly  from  spectral  information. 
MSAVI  also  increases  the  sensitivity  to  vegetation  and  minimizes  soil  influences 
(Qi  et  al.  1994).  The  MSAVT  accounts  for  possible  variations  of  soil  reflectance 
through  an  inductive  method  based  on  the  opposite  trends  of  NDVI  and  WDVI, 
thus  eliminating  subjective  soil  correction  L  assignments.  In  the  original  field 
test  in  cotton  fields,  density  of  the  cotton  canopy  was  predicted  more  accurately 
by  MSAVI  than  other  derivations  of  SAVE  using  remotely  sensed  imagery.  Not 
only  does  MSAVI  eliminate  the  requirement  to  estimate  L,  but  for  arid  lands, 
MSAVI  has  been  promoted  as  the  best  predictor  of  shrub  cover  measurements 
among  those  indices  that  attempt  to  mitigate  the  influence  of  soil  reflectance 
(Rondeaux  and  Steven  1996;  Senseman,  Bagley,  and  Tweddale  1996). 


Brightness  Indices 

Albedo  (or  reflectance)  and  soil  brightness  indices  have  also  been  calculated  for 
the  purpose  of  estimating  percent  vegetative  cover  from  satellite  imagery.  While 
vegetation  or  greenness  indices  attempt  to  measure  the  reflectance  of  vegetation 
directly  while  minimizing  soil  effects,  brightness  indices  measure  the  total 
brightness  or  reflectance  of  the  Earth’s  surface,  and  therefore  are  more  sensitive 
to  soil  background  reflectance.  Vegetation  tends  to  mask  the  reflectance  of  soils, 
particularly  in  sparsely  vegetated  arid  environments  with  highly  reflective  back¬ 
ground  soils.  Therefore,  an  inverse  relationship  between  brightness  indices  and 
vegetative  cover  has  been  used  to  estimate  vegetative  cover. 

Albedo  is  defined  as  the  ratio  of  all  shortwave  radiation  reflected  by  the  Earth's 
surface  to  solar  irradiance  incident  on  the  surface.  Albedo  has  been  correlated 
with  measurements  of  vegetative  cover  in  several  arid  and  semi-arid  environ¬ 
ments  (Robinove  et  al.  1981;  Musick  1986;  Frank  1984).  Planetary  Albedo,  or 
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planetary  reflectance,  is  a  measure  of  reflectance  of  the  entire  Earth-atmosphere 
system  and  is  calculated  directly  from  observations  recorded  at  the  satellite. 
Surface  reflectance,  or  surface  Albedo,  is  a  measure  of  reflectance  of  the  Earth’s 
surface.  To  calculate  surface  Albedo  or  reflectance  using  remotely  sensed  obser¬ 
vations,  radiometric  and  atmospheric  corrections  must  be  applied  to  correct  for 
atmosphere,  topography,  and  solar  geometry  variations.  Many  times,  reflectance 
or  Albedo  is  calculated  across  all  of  the  visible  and  near  infrared  wavelengths 
available.  Other  times,  reflectance  of  a  single  spectral  band,  or  in-band  Albedo, 
often  correlates  well  with  vegetative  cover  in  arid  environments.  In  all  cases, 
vegetative  cover  is  measured  indirectly  by  the  relative  decrease  in  brightness  of 
the  background  soils  due  to  vegetative  cover  or  shadowing. 

Other  soil  brightness  indices  have  also  been  developed  to  assess  soil  brightness, 
and  indirectly,  vegetation  cover.  One  common  index  is  the  Kauth-Thomas  Soil 
Brightness  Index  (SBI,  Kauth  et  al.  1978).  The  SBI  is  calculated  based  on  the 
Tasseled  Cap  transformation  described  earlier.  Similar  to  the  Greenness  Index, 
the  SBI  is  derived  from  the  characteristic  plot  of  red  vs.  near  infrared  reflectance 
in  spectral  space.  Soil  brightness  indices  have  been  correlated  with  arid  shrub 
cover  in  several  studies  at  the  Jornada  Long  Term  Ecological  Research  project, 
which  is  within  close  proximity  to  Fort  Bliss  (Duncan  et  al.  1993;  Mustek  1984, 
1986). 

Change  Detection 

Remotely  sensed  imagery  provides  an  ideal  data  source  for  assessing  change  in 
surface  condition  because  it  has  high  temporal  frequency  and  large  aerial  cover¬ 
age.  Several  techniques  have  been  developed  to  assess  change  using  temporal 
image  data  sets  (Singh  1989).  Some  techniques  may  work  better  than  others  in 
different  ecological  settings.  In  arid  and  semi-arid  environments  such  as  Fort 
Bliss,  the  most  important  indicator  of  land  condition  is  vegetative  cover.  Vegeta¬ 
tion  Indices,  or  Greenness  Indices,  and  to  a  greater  extent.  Brightness  Indices, 
have  exhibited  high  correlation  with  vegetative  cover  in  arid  environments. 
Typically,  relative  increases  in  soil  brightness  or  Albedo  are  generally  associated 
with  degradation,  while  relative  decreases  are  generally  associated  ivith  im¬ 
provement  of  land  condition.  Conversely,  relative  increases  in  vegetation  indices 
are  associated  with  increased  vegetative  cover,  and  therefore  improved  land  con¬ 
dition,  while  decreases  in  vegetation  indices  indicate  a  loss  of  vegetative  cover, 
and  therefore,  high  erosion  potential.  In  many  cases,  a  threshold  is  established 
that  represents  “normal”  change  in  indices,  and  areas  that  increase  or  decrease 
beyond  this  threshold  are  identified  as  improved  or  degraded  (Musick  1984, 
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1986;  Robinove  et  al.  1981;  Frank  1984;  Price,  Pyke,  and  Mendes  1992;  and  Yool, 
Makaio,  and  Watts  1997). 

Simple  observations  of  trends  or  differencing  of  indices  from  different  image 
dates  provides  a  relative  measure  of  change  in  cover.  Other  variables  may  affect 
trends  in  indices,  such  as  soil  moisture,  climate,  and  phenology  of  vegetation.  In 
addition,  technical  errors  in  calibration  and  georegistration  of  temporal  imagery 
can  also  bias  results.  However,  because  complete  field  surveys  are  impossible 
due  to  cost  and  time  constraints,  monitoring  change  in  indices  still  provides  a 
cost-effective  method  to  monitor  general  and  relative  trends  in  land  condition  for 
military  installations. 
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3  Study  Site 

A  rectangular  study  site  approximately  21  km  by  9  km  was  chosen  in  the  Tu- 
larosa  Valley  at  Fort  Bliss  (Figure  1).  U.S.  Route  54  divides  the  study  site  into 
East  and  West  subsections.  The  western  edge  of  the  study  site  begins  near  Scott 
Tank  and  extends  21  km  east,  crossing  U.S.  Route  54  and  extending  into  the 
McGregor  Guided  Missile  Range.  The  site  was  chosen  to  represent  a  contrast  of 
historic  training  impact.  The  section  west  of  U.S.  Route  54  (hereafter  called 
West)  is  in  the  Fort  Bliss  Dona  Ana-Orogrande  Complex  and  includes  Maneuver 
Areas  5C,  5D,  and  parts  of  4B  and  4C.  This  area  has  been  heavily  impacted  by 
military  training  and  maneuvering,  with  very  sparse  vegetation  remaining  in 
the  interdunal  spaces  between  the  coppice  dunes.  The  area  to  the  east  of  U.S. 
Route  54  (hereafter  called  East)  is  in  the  McGregor  Guided  Missile  Range.  This 
section  has  not  been  impacted  by  tracked  and  wheeled  vehicle  msmeuvers,  and 
therefore  serves  as  a  control  for  comparison  with  the  heavily  impacted  maneuver 
areas. 

The  average  elevation  of  Tularosa  Valley  in  and  around  the  study  site  is  ap¬ 
proximately  4000  ft  (1220  meters  [m]).  The  study  area  is  comprised  primarily  of 
mesquite-covered  coppice  dunes.  The  dunes  are  typically  oblong  and  average 
11  m  long  and  6  m  wide,  with  average  interdunal  spacing  of  8  m.  Dunes  are  pre¬ 
dominantly  covered  by  mesquite  and  occasional  fourwing  saltbush,  ranging  in 
height  from  1  to  4  m  with  an  average  height  of  2  m  and  crown  density  from  35  to 
100  percent,  averaging  55  percent.  Interdunal  vegetation  is  dominated  by  broom 
snakeweed,  plains  lovegrass,  fluffgrass,  whitestem  paperflower,  and  isolated 
creosote  bush  (Fort  Bliss  Terrain  Analysis  [map  book]  1978). 
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4  Data  Processing  and  Methodology 

Satellite  Imagery 

Landsat  Multi-Spectral  Scanner  (MSS)  and  Landsat  Thematic  Mapper  (TM)  im¬ 
agery  was  compiled  and  analyzed  for  this  study.  The  MSS  sensor  was  launched 
aboard  the  first  civil  remote  sensing  satellite,  now  known  as  Landsat  1,  in  1972. 
MSS  collected  reflectance  in  four  (later  five)  spectral  bands  ranging  from  blue  to 
near-infrared  (0.5  to  1.1  |im)  at  an  average  spatial  resolution  of  80  m  (Table  1). 
Two  more  MSS  sensors  were  launched  aboard  Landsat  2  in  1975  and  Landsat  3 
in  1978,  each  with  improved  sensor  and  communication  capabilities.  In  1982,  the 
TM  sensor  was  launched  along  with  MSS  on  Landsat  4.  TM  recorded  reflected 
and  emitted  energy  in  seven  narrower  spectral  bands  from  blue  to  thermal  infra¬ 
red  (0.45  jxm  to  11.66  |J.m).  A  second  TM  sensor  was  launched  aboard  Landsat  5 
in  1984.  TM  sensors  are  still  in  operation  and  images  are  collected  daily  arotmd 
the  world. 


Table  1.  Wavelength  characteristics  of  Landsat  MSS  and  Landsat  TM  sensors*. 


Channel 

Band  width  (micrometers) 

Ground  IFOV**  (m) 

Landsat  1  through  5  MSS 

1 

0.5  -  0.6 

80 

2 

0.6  -  0.7 

80 

3 

0.7  -  0.8 

80 

4 

0.8-  1.1 

80 

Landsat  4  and  5  TM 

1 

0.45  -  0.52 

30 

2 

0.53  -  0.60 

30 

3 

0.63  -  0.69 

30 

4 

0.76  -  0.90 

30 

5 

1.55-  1.75 

30 

6 

10.42-12.50 

120 

7 

2.08  -  2.35 

30 

*  Adopted  from  Table  9-2  of  Lillesand  and  Kiefer  1987. 
**  Instantaneous  field  of  view. 
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Sensor  Differences 

Although  both  MSS  and  TM  sensors  collect  reflectance  in  the  red  and  near- 
infrared  regions  of  the  electromagnetic  spectrum,  there  is  not  a  direct  correspon¬ 
dence  between  the  bandwidths  that  are  collected  by  each  sensor.  MSS  Band  2  is 
similar  in  spectral  coverage  to  TM  Band  3,  and  these  two  bands  were  used  as  in¬ 
put  for  red  reflectance  in  greenness  and  brightness  indices.  MSS  Band  4  is 
roughly  similar  to  TM  Band  4,  although  the  MSS  bandwidth  is  much  greater. 
Band  4  from  both  sensors  was  used  as  input  for  near-infrared  reflectance  in 
greenness  and  brightness  indices. 


Image  Compilation 

Imagery  was  acquired  from  several  sources  and  different  sensors  to  compile  an 
imagery  data  set  covering  a  23-year  time  period  from  1972  through  1994.  An 
initial  review  of  the  existing  Fort  Bliss  imagery  already  available  at  the  Engi¬ 
neer  Research  and  Development  Center/Construction  Engineering  Research 
Laboratory  (ERDC/CERL)  and  the  University  of  Wisconsin-Madison  was  docu¬ 
mented.  Three  images  of  Fort  Bliss  were  previously  purchased  by  the  UW- 
Madison  1994-1995  “Environmental  Monitoring  Practicum,”  which  was  funded 
by  the  U.S.  Army  Environmental  Center.  Temporal  gaps  in  coverage  of  existing 
imagery  were  identified  and  an  archival  search  of  available  imagery  for  these 
time  periods  was  completed.  Additional  imagery  was  acquired  to  provide  a  sin¬ 
gle  continuous  temporal  data  set  of  19  images  spanning  from  1972  through  1994. 
The  late  summer  and  early  autumn  (August/September/October)  was  identified 
as  a  target  anniversary  date  for  temporal  analysis.  Peak  precipitation  generally 
occurs  during  this  time  period,  which  in  ttu*n  results  in  peak  vegetative  growth, 
vigor,  and  photosynthetic  activity.  However,  if  there  was  existing  imagery  for 
Fort  Bliss,  it  was  included  for  analysis,  even  if  the  acquisition  date  did  not  match 
this  target  anniversary  date.  As  a  result,  some  image  acquisition  dates  are  dis¬ 
persed  around  this  target  anniversary  date,  and  inclusion  of  dates  that  vary  sig¬ 
nificantly  from  the  target  anniversary  date  are  noted  as  such  and  are  not  in¬ 
cluded  in  analysis  of  temporal  trends  in  land  condition.  To  conduct  a  cross 
comparison  of  TM  and  MSS  sensors,  separate  TM  and  MSS  scenes  of  Fort  Bliss, 
both  acquired  on  28  October  1984,  were  included  in  the  study.  Table  2  lists  the 
dates  and  scene  identification  numbers  of  MSS  and  TM  images  used  in  this 
study. 
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Table  2.  MSS  and  TM  images  used  in  this  study. 


MSS 

Date 

Scene  ID 

6  September  1972 

01196103100700001 

23  December  1972 

LM1 035038007235890 

20  October  1 974 

LM1 035038007429390 

25  November  1974 

LM1 035038007432990 

30  September  1976 

LM2035038007627490 

29  September  1 978 

LM3035038007827290 

26  December  1 980 

LM2035038008036190 

28  June  1982 

01196103100700002 

5  December  1983 

LM4033038008333990 

28  October  1 984 

LM5033038008430290 

T 

M 

28  October  1984 

LT5033038008430210 

1 5  August  1 986 

LT5033038008622710 

19  September  1987 

LT5033038008726210 

24  October  1 988 

LT5033038008829710 

29  October  1 990 

LT5033038009028610 

4  June  1992 

LT4033038009215610 

1 8  October  1 982 

LT5033038009229210 

19  September  1993 

LT5033038009326210 

9  November  1 994 

LT5033038009431310 

Calibration 

Reflectance  recorded  by  satellites  is  influenced  not  only  by  the  spectral  proper¬ 
ties  of  the  Earth’s  surface,  but  also  by  the  Earth’s  atmosphere,  solar  geometry, 
topographic  relief,  and  instrument  design  and  deterioration.  lb  compare  images 
collected  at  different  times  or  by  different  sensors,  these  effects  must  be  mini¬ 
mized.  The  process  by  which  dimensionless  digital  numbers  are  converted  to 
physical  units  of  reflectance  is  called  radiometric  correction  or  calibration.  Sev¬ 
eral  radiometric  calibration  techniques  of  varying  complexity  can  be  applied.  In 
many  cases,  complete  radiometric  calibration  requires  field  measurements  of 
atmospheric  influences  at  the  time  of  image  acquisition.  In  addition,  information 
on  solar  zenith  angle,  topography  of  the  Earth’s  surface,  and  radiometric  quali¬ 
ties  of  the  sensor  are  needed.  Solar  zenith  angle  and  radiometric  information 
are  recorded  in  the  image  header  files.  Prelaunch  radiometric  sensitivity  of  the 
sensor  is  also  published  information.  For  this  study,  in-situ  measurements  of 
atmospheric  scattering  were  unavailable.  Therefore,  an  alternative  method  was 
used  to  estimate  and  correct  for  atmospheric  scattering  in  the  temporal  data  set. 
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A  large  pond  in  Sunlake  Park  near  El  Paso,  Texas,  was  selected  to  represent  a 
relatively  stable  target  over  the  entire  period  from  1972  through  1994.  For  each 
image,  the  mean  value  for  all  the  pixels  within  the  area  defined  by  Sunlake  Park 
pond  were  calculated  for  each  band  (see  Appendix  A).  In  the  MSS  images,  the 
means  were  relatively  stable  for  all  bands  for  all  years.  For  the  TM  images,  the 
means  were  relatively  stable  in  bands  1  through  4  but  varied  greatly  for  bands  5 
and  6. 

Any  differences  in  observed  reflectance  values  for  the  pond  from  different  dates 
were  assumed  to  be  due  to  a  change  in  atmospheric  attenuation.  Using  these 
mean  values  for  the  pond  for  each  band,  the  minimum  mean  for  each  band  was 
identified  and  used  as  a  standard  or  baseline.  This  minimum  mean  was  sub¬ 
tracted  from  the  mean  values  observed  in  the  same  band  for  the  other  dates. 
The  resulting  differences,  or  offsets,  were  assumed  to  be  the  atmospheric  influ¬ 
ence  present  in  each  image.  These  offsets  for  each  band  were  then  subtracted 
from  all  pixel  values  in  the  respective  images,  resulting  in  an  atmospherically 
corrected  temporal  data  set  (Scarpace  et  al.  1995).  It  is  recogmzed  that  a  natural 
water  body  is  not  an  ideal  stable  target  for  the  purpose  of  image  calibration,  as 
there  may  have  been  variations  in  reflectance  of  the  pond  due  to  depth  and  tur¬ 
bidity  of  the  water.  However,  there  were  limited  stable  surfaces  that  could  be 
identified  in  each  scene  of  the  temporal  data  set.  Therefore,  the  pond  was  chosen 
for  calibration  purposes,  and  any  variations  in  pond  reflectance  due  to  actual 
changes  in  the  lake  itself  were  assumed  to  be  minimal. 

Solar  illumination  differences  due  to  aspect  from  varying  topography  were  as¬ 
sumed  to  be  minimal  for  the  study  site,  and  therefore  were  not  considered  in  the 
analysis.  Solar  zenith  angles  were  acqtured  firom  image  header  files.  When 
available,  radiometric  offset  and  gain  values  were  also  acquired  from  image 
headers.  If  unavailable,  published,  pre-launch  radiometric  offset  and  gain  val¬ 
ues  were  used  to  calculate  reflectance. 


Noise  Reduction 

Systematic  horizontal  strips  of  bad  data  were  evident  in  several  of  the  images. 
These  problems  were  particularly  evident  in  the  MSS  imagery,  although  some 
noise  is  also  apparent  in  the  TM  imagery.  Striping  of  MSS  imagery  is  a  result  of 
a  known  problem  with  improper  synchronization  of  scan  lines  by  the  sensor.  A 
periodic  noise  removal  algorithm  (available  in  Earth  Resources  Data  Analysis 
Systems,  Inc.  [ERDAS])  Imagine  software)  was  applied  to  striped  images  in  an 
attempt  to  remove  or  minimize  noise.  The  periodic  noise  removal  performs  an 
enhancement  of  the  image  based  on  a  Fourier  transform.  The  original  striped 
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image  is  first  divided  into  128  by  128  pixel  blocks.  The  Fourier  transform  for 
each  block  is  calculated  and  the  log-magnitudes  of  each  Fast  Fourier  Transform 
(FFT)  block  are  averaged.  The  averaging  removes  or  significantly  reduces  all  of 
the  periodic  interference  (ERDAS  1993).  However,  the  periodic  noise  removal 
changed  the  actual  digital  numbers  in  the  imagery,  and  therefore  had  a  minor 
effect  on  vegetation  and  soil  brightness  indices. 


Georeferencing 

Once  imagery  was  compiled  and  corrected  for  atmospheric  and  problem  scan 
lines,  the  images  were  georeferenced  to  a  Universal  Transverse  Mercator  (UTM) 
coordinate  system.  Satellite  imagery  is  acquired  in  raw  form,  and  therefore  is 
not  referenced  to  an  established  coordinate  system.  Georectification  is  required 
to  locate  reflectance  values  of  individual  pixels  to  a  known  coordinate  system. 
Once  an  image  is  georeferenced,  reflectance  values  can  be  related  to  a  geographic 
location  on  the  ground.  Typically,  georeferenced  imagery  is  entered  as  a  data 
layer  in  a  geographic  information  system  (GIS)  for  subsequent  spatial  analysis 
with  other  geographic  information,  including  additional  imagery.  To  georefer¬ 
ence  an  image,  precise  geographic  coordinates  of  identifiable  ground  control 
points  (GCPs)  must  be  located  in  the  imagery.  Geographic  coordinates  for  con¬ 
trol  points  can  be  obtained  by  locating  the  same  points  in  another  map  that  is 
already  georeferenced  to  a  coordinate  system,  or  the  coordinates  can  be  collected 
in  the  field  using  Global  Positioning  System  (GPS)  technology.  Based  on  known 
control  points,  a  transformation  matrix  is  calculated  and  used  to  resample  each 
pixel  or  data  element  to  a  known  coordinate  system. 

The  University  of  Wisconsin-Madison  1994-1995  “Environmental  Monitoring 
Practicum”  collected  coordinates  for  ground  control  points  using  GPS  technology 
and  georeferenced  the  4  June  1992  image  to  (1)  a  UTM  coordinate  system  and 
projection,  (2)  Zone  13,  and  (3)  Northern  Hemisphere  Clarke  1866  Spheroid.  All 
other  images  were  georeferenced  to  this  base  image.  Although  some  images  of 
Fort  Bliss  had  been  georeferenced  prior  to  this  study,  the  registration  was  of 
questionable  quality.  Therefore,  original  source  imagery  was  located,  control 
points  were  identified,  and  the  images  were  georeferenced  to  the  base  1992  TM 
image. 

Stable  GCPs  such  as  road  intersections  and  buildings  were  located  in  the  raw 
source  image  and  the  base  1992  image.  Once  a  suitable  number  of  GCPs  were 
selected,  a  transformation  matrix  and  corresponding  measure  of  fit,  called  the 
Root  Mean  Square  (RMS)  error,  was  calculated.  A  transformation  matrix  is  a  set 
of  coefficients  used  to  convert  pixels  to  georeferenced  coordinates.  RMS  error  is 
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the  distance  between  the  desired  output  coordinate  for  a  GCP  and  the  actual 
output  coordinate  for  the  same  point  if  that  point  were  to  be  transformed  with 
the  current  transformation  matrix.  Therefore,  when  selecting  and  evaluating 
GCPs,  the  goal  was  to  minimize  RMS  error  while  distributing  GCPs  evenly 
across  the  image.  RMS  error  can  be  measured  in  meters  on  the  ground  or  in  a 
fraction  of  the  smallest  picture  element  (pixel)  in  the  image.  Using  TM  imagery 
with  30-m  spatial  resolution  as  an  example,  an  RMS  error  of  .333  would  indicate 
that  the  base  and  source  images  are  spatially  misregistered  by  one-third  of  a 
pixel  (or  10  meters).  A  maximum  acceptable  RMS  error  of  .30  was  established 
for  all  georectifications  in  this  study. 

An  attempt  was  made  to  use  the  same  GCPs  in  successive  images  when  possible. 
However,  over  the  time  period  of  this  study,  many  landmarks  that  were  clearly 
identifiable  in  the  1992  base  image  could  not  be  located  in  imagery  from  the 
early  1970s.  Therefore,  new  GCPs  had  to  be  located;  the  same  control  points 
were  not  used  for  all  images  in  the  temporal  data  set. 

Once  an  acceptable  set  of  GCPs  was  located  in  the  imagery,  a  first  order  polyno¬ 
mial  transformation  was  applied  and  the  raw  image  pixels  were  resampled  to 
match  the  1992  image  using  a  nearest  neighbor  resampling  in  ERDAS  Imagine. 

Final  evaluation  involved  visual  comparison  of  the  newly  georeferenced  image 
with  the  1992  base  image.  In  the  event  that  georegistration  appeared  inaccu¬ 
rate,  GCPs  with  unusually  high  RMS  errors  were  eliminated  and  different  con¬ 
trol  points  were  selected.  This  process  was  repeated  until  low  RMS  errors  were 
reported  and  visual  inspection  confirmed  an  accurate  georectification. 


Study  Site  Subset 

Once  all  images  in  the  temporal  data  set  were  georeferenced  to  a  common  coor¬ 
dinate  system,  an  area  corresponding  to  the  study  site  was  extracted  from  each 
image.  For  MSS  imagery,  the  resulting  image  subsets  matching  the  study  area 
were  259  by  109  pixels,  while  subsets  from  the  TM  imagery  were  691  by  291  pix¬ 
els. 


Index  Calculations 

Four  vegetation  or  greenness  indices  and  four  brightness  indices  were  calculated 
for  each  of  19  images  in  the  temporal  data  set,  which  spanned  from  1972  through 
1994.  Each  index  was  selected  based  on  its  sensitivity  to  percent  vegetative 
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cover  as  summarized  in  previous  literatrure  (Tucker  1979;  Robinove  et  al.  1981; 
McDaniel  and  Haas,  1982;  Musick  1984;  Huete  1988;  Duncan  et  al.  1993;  Qi  et 
al.  1994).  The  vegetation  or  greenness  indices  calculated  were  the  Normalized 
Difference  Vegetation  Index  (NDVI),  the  “Tasseled  Cap”  Greenness  Index,  and 
two  derivations  of  the  Soil-Adjusted  Vegetation  Index:  SAVI  and  a  Modified 
SAVI  (MSAVI2).  The  brightness  indices  calculated  were  Albedo  (or  reflectance) 
over  all  visible  and  near-infrared  bands,  in-band  Albedo  in  the  red  and  near- 
infrared  wavelengths  individually,  and  the  “Tasseled  Cap”  Soil  Brightness  Index. 
All  processing  was  done  under  the  image  processing  software  package  Imagine, 
which  is  distributed  by  ERDAS.  A  graphical  model  in  Imagine  was  used  to  cal¬ 
culate  the  indices  as  follows: 

1.  NDVI  (Rouse  et  al.  1974) 


NDVI  = 


NIR-red 
NIR  +red 


2.  SAVI  (Huete  1988) 


SAVI  = 


NIR-red 
NIR  -\-red  -\-L 


(1+Z.) 


3.  MSAVI2  (Qi  et  al.  1994) 


MSAVI2  =-(2(NIR  +1) -^(2 xNIR  +iy  -8( M/?  -red )) 


4.  Tasseled  Cap  Greenness  (TCG)  Vegetation  Index 
For  MSS: 

TCG  =  (-0.283)(MSS1)  -  (0.660)(MSS2)  +  (0.577)(MSS3)  +  (0.388)(MSS4) 

For  TM4: 

TCG  =  (■0.285)(TM1)  ■  (.244)(TM2)  ■  (.544)(TM3)  +  (.724)(TM4)  +  (.084)(TM5)  -  (0.180)(TM7) 
For  TM5: 


TCG  =  (-0.273)(TM1)  -  (.217)(TM2)  -  (.551)(TM3)  +  (.722)(TM4)  +  (.073)(TM5)  -  (0.165)(TM7) 
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5.  Albedo 

Prior  to  calculation  of  Albedo  (or  reflectance),  digital  numbers  were  converted  to 
spectral  band  radiance  for  MSS  and  TM  imagery  using  (Markham  and  Barker 
1986): 

Lmm  A+[Lmax  A- Lmin  x^N 

D  max  Equation  1 

where:  Lmin  =  spectral  radiance  of  each  band  at  DN  =  0  in 

mWcm  ^  Sr  ^:m  ^ 

Lmax  =  spectral  radiance  of  each  band  at  DN=255  (DN=  127  for  MSS)  in 
mWcm^Sr  '^:m  ^ 


Dmax  =  range  of  rescaled  radiance  in  DN 
DN  =  input  digital  number. 

Once  spectral  radiance  was  calculated  according  to  Equation  1,  Albedo  was  cal¬ 
culated  as  the  ratio  of  reflected  solar  radiation  to  incoming  solar  irradiance: 

^~~E^  Equation  2 

-2  -1  -1 

where:  Lg  =  spectral  radiance  in  mWcm  Sr  :m  (from  Equation  1) 
d  =  Earth  -  sun  distance  in  astronomical  units 

Eg  =  exoatmospheric  spectral  irradiance  at  the  top  of  the  atmosphere 
in  mWcm"^  :m  ^ 

2  =  solar  zenith  angle. 

Radiance  is  derived  from  digital  numbers  recorded  on  computer-compatible 
tapes.  Exoatmospheric  spectral  irradiances  at  the  top  of  the  atmosphere  are  es¬ 
timates  derived  from  Markham  and  Barker  (1986).  Correcting  exoatmospheric 
solar  irradiance  in  the  denominator  by  cos  2  normalizes  scenes  to  an  overhead  or 
nadir  sun  angle  and  accounts  for  differences  in  solar  irradiance  for  the  time  of 
day  and  day  of  the  year  (Robinove  et  al.  1981). 
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Total  reflectance  or  planetary  Albedo  was  calculated  by  integrating  spectral  ra¬ 
diance  and  irradiances  across  all  of  the  visible  and  the  near-  and  middle-infrared 
bands.  Planetary  Albedo  for  MSS  was  integrated  across  all  4  bands,  while 
planetary  Albedo  for  TM  imagery  was  integrated  across  only  bands  1  through  5. 
In-band  planetary  Albedo  for  red  and  near-infrared  bands  was  calculated  in  the 
same  manner  as  total  reflectance  or  planetary  Albedo  (Equation  1).  However, 
red  and  near-infrared  reflectances  were  calculated  using  radiance  and  irradi¬ 
ances  in  only  the  red  and  near-infrared  wavelengths. 

6.  Tasseled  Cap  Brightness  (TCB)  Index 

For  MSS: 

TCB  =  (.332)(MSS4)  +  (.603)(MSS5)  +  (.675)(MSS6)  +  (.262)(MSS7) 

For  TM4: 

TCB  =  (.304)(TM1)  +  (.279)(TM2)  +  (.474)(TM3)  +  (.559)(TM4)  +  (.508)(TM5)  +  (.186)  (TM7) 
For  TM5: 


TCB  =  (.291)(TM1)  +  (.249)(TM2)  +  (.481)(TM3)  +  (.557)(TM4)  +  (.444)(TM5)  +  (.171)  (TM7) 
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5  Analysis  and  Results 

A  temporal  data  set  of  19  satellite  images  spanning  23  years  was  compiled  and 
analyzed  to  assess  historical  trends  in  relative  land  condition  and  contrasting 
land  uses  over  this  time  period  for  a  specific  study  site  at  Fort  Bliss,  Texas. 
However,  a  secondary  product  of  this  research  is  a  complete  temporal  data  set  for 
the  entire  installation  over  this  same  time  period.  Spanning  from  1972  through 
1994,  the  data  set  includes  imagery  ranging  from  the  earliest  commercially 
available  satellite  platform  (Landsat  1  MSS)  to  a  satellite  platform  (Landsat  5 
TM)  that  is  still  functional  today.  For  most  years  in  this  temporal  data  set,  com¬ 
plete  coverage  of  the  entire  installation  does  exist.  Many  images  were  compiled 
from  archival  information  at  ERDC/CERL,  while  others  were  purchased  as  part 
of  this  project.  Because  of  the  size  and  location  of  Fort  Bliss,  several  images  are 
required  to  provide  complete  geographic  coverage  of  the  installation.  For  exam¬ 
ple,  three  Landsat  TM  scenes  are  required  for  complete  coverage  because  of  the 
orbit  of  the  Landsat  satellite  and  swath  size  of  the  TM  sensor.  However,  two  TM 
scenes  provide  close  to  95  percent  coverage  of  the  installation.  A  third  image  is 
required  to  cover  the  extreme  northeast  comer  of  Fort  Bliss,  and  in  many  cases, 
was  not  available  from  archival  imagery  of  the  installation.  For  the  years  in 
which  imagery  was  acquired  specifically  for  this  project,  a  single  TM  scene  was 
selected,  which  provided  coverage  of  the  study  site.  Therefore,  for  those  years, 
complete  coverage  of  the  entire  installation  is  not  available. 

The  results  summarized  in  this  chapter  are  specifically  focused  on  a  retrospec¬ 
tive  analysis  of  a  specific  study  site,  as  described  in  Chapter  3,  and  were  not 
compiled  for  the  entire  installation.  The  analysis  is  focused  on  historical  trends 
and  spatial  patterns  of  remotely  sensed  indices,  which  are  indicators  of  vegeta¬ 
tive  cover,  and  indirectly,  land  condition  for  the  entire  study  site,  and  also  fo¬ 
cused  on  areas  of  contrasting  land  use  within  the  study  site. 


Historical  Trends 

Eight  remotely  sensed  indices  were  calculated  from  each  image  date  for  the 
study  site.  Four  of  the  indices  were  vegetation  indices  or  greenness  indices; 
(1)  the  Normalized  Difference  Vegetation  Index  (NDVI),  (2)  the  “Tasseled  Cap” 
Greenness  Index  (TCG),  (3)  the  Soil-Adjusted  Vegetation  Index  (SAVI),  and 
(4)  the  Modified  SAVE  (MSAVI2).  The  four  remaining  indices  were  brightness 
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indices:  (1)  Albedo,  (2)  the  “Tasseled  Cap”  Brightness  Index  (TCB),  (3)  Red  Re¬ 
flectance,  and  (4)  Near-Infrared  Reflectance. 

When  analyzing  a  single  image,  radiometric  and  atmospheric  calibration  was  not 
necessary.  However,  because  the  primary  goal  was  to  analyze  a  temporal  se¬ 
quence  of  images  collected  from  different  platforms  and  under  presumably  dif¬ 
ferent  atmospheric  conditions,  radiometric  and  atmospheric  corrections  were  ap¬ 
plied  to  the  original  Digital  Numbers  (DNs),  and  all  indices  were  calculated 
using  calibrated  reflectance  data.  The  resulting  indices  calculated  from  cali¬ 
brated  reflectance  values  were  dimensionless  numbers.  Therefore,  the  indices 
were  ideal  for  comparisons  between  pixels  in  a  single  image,  or  in  this  case, 
when  calibrated,  for  comparison  between  pixels  representing  the  same  geo¬ 
graphic  area  in  a  temporal  sequence.  Similarly,  reflectance  in  the  red  and  near- 
infrared  wavelengths  was  also  calibrated  for  comparison  over  the  same  temporal 
sequence.  However,  radiometric  and  atmospheric  calibration  is  not  an  exact  sci¬ 
ence,  and  although  procedures  were  followed  to  calibrate  image  data  from  differ¬ 
ent  dates  for  the  purpose  of  isolating  actual  differences  of  the  Earth’s  surface, 
the  results  should  be  used  only  for  identif5dng  general  and  relative  trends  in  sur¬ 
face  conditions. 

Especially  problematic  to  this  research  was  the  fact  that  imagery  was  compiled 
from  Landsat  MSS  and  TM  sensors,  which  have  significantly  different  spatial 
and  spectral  characteristics.  Due  to  differences  between  these  sensors,  compari¬ 
son  of  index  values  spanning  between  MSS  and  TM  dates  may  not  provide  an 
accurate  representation  of  actual  surface  conditions  over  this  entire  time  period. 
Because  different  spatial  and  spectral  resolution  imagery  captures  a  different 
mixture  of  landcover  components  within  a  single  pixel  (or  data  element),  the  re¬ 
sulting  indices  calculated  from  MSS  data  are  not  suitable  for  comparison  with 
indices  calculated  from  TM  data.  Therefore,  the  temporal  sequence  of  MSS  im¬ 
agery  spanning  from  1972  through  1984  was  analyzed  and  is  presented  sepa¬ 
rately  from  the  temporal  sequence  of  TM  imagery,  which  spans  from  1984 
through  1994. 

Appendix  B  lists  the  mean,  standard  deviation,  minimum  and  maximum  index 
values,  and  reflectance  values  for  the  entire  study  site  for  each  image  date  and 
for  each  of  the  eight  greenness  and  brightness  indices.  Appendix  B  also  includes 
similar  statistics  for  subsets  of  the  study  area,  which  are  discussed  later  in  this 
report.  Although  eight  indices  were  tabulated  for  each  date,  a  final  group  of  four 
indices  (Albedo,  Tasseled  Cap  Brightness,  Red  Reflectance,  and  MSAVI2)  were 
selected  for  more  detailed  analysis.  These  indices  were  selected  based  on  related 
research  and  literature  (Musick  1984;  Robinove  et  al.  1981;  Duncan  et  al.  1993; 
Pickup,  Chewings,  and  Nelson  1993),  which  indicated  that  they  would  be  most 
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sensitive  to  vegetative  cover  and  land  condition,  and  therefore,  would  be  optimal 
for  analyzing  trends  in  land  condition  over  time.  Mean  index  values  for  the  en¬ 
tire  study  site  for  each  of  these  four  indices  were  plotted  for  each  image  date  in 
Figures  2  and  3. 

Initial  examination  of  the  three  mean  brightness  indices  (Albedo,  Red  Reflec¬ 
tance,  and  “Tasseled  Cap”)  from  1984  through  1994  using  TM  imagery  (Figure  3) 
reveals  a  similar  pattern.  Each  index  dropped  sharply  in  both  1986  and  1992. 
However,  note  that  the  1992  image  was  collected  in  June,  which  varies  consid¬ 
erably  from  the  remaining  acquisition  dates  that  span  from  August  through  Oc¬ 
tober.  Ideally,  when  evaluating  a  temporal  sequence  of  imagery  to  detect  change, 
images  should  be  acquired  as  close  as  possible  to  the  same  anniversary  date  to 
minimize  phenological  differences  between  scenes.  Therefore,  when  evaluating 
trends  over  this  time  period,  the  June  1992  data  should  be  interpreted  with  cau¬ 
tion.  The  1986  scene  was  collected  in  August,  which  represents  the  earliest  part 
of  the  targeted  image  acquisition  time  period.  Therefore,  phenological  differences 
may  also  explain  the  relatively  low  brightness  indices  for  this  date  (Albedo  = 
11.91,  TCB  =  44.64,  Red  =  11.89),  yet  the  indices  are  still  considerably  lower  than 
the  1988  scene  (Albedo  =  14.22,  TCB  =  52.52,  Red  =  14.45),  which  was  also  ac¬ 
quired  in  August. 

Another  critical  consideration  when  evaluating  trends  in  spectral  indices  to  as¬ 
sess  land  condition  is  climatic  information,  such  as  antecedent  precipitation  lead¬ 
ing  up  to  the  image  acquisition  date.  Precipitation  affects  spectral  brightness 
indices  in  two  ways.  First,  any  precipitation  that  occurs  close  to  the  acquisition 
date  of  the  image  may  result  in  increased  soil  moisture  that  may  act  to  decrease 
brightness,  as  wet  or  damp  soils  are  darker  than  dry  soils.  Second,  increased 
precipitation  usually  results  in  an  increase  in  the  amount  and  vigor  of  perennials 
and  annuals,  which  may  also  result  in  decreased  brightness  (Robinove  et  al. 
1981).  Robinove  and  associates  revealed  that  the  major  determining  factor  in  a 
change  in  brightness  measures  such  as  Albedo,  for  example,  was  antecedent  pre¬ 
cipitation  for  7  or  more  months  leading  up  to  the  later  date  in  the  temporal  se¬ 
quence  over  which  change  in  Albedo  is  observed.  These  results  have  not  been 
validated  with  rigorous  statistical  analysis,  but  they  do  demonstrate  the  impor¬ 
tance  of  considering  climatic  variables  when  analyzing  change  in  spectral  indi¬ 
ces.  Therefore,  monthly  precipitation  data  was  collected  from  1972  through  1995 
from  El  Paso  Airport,  which  is  a  few  kilometers  south  of  the  study  area  (Figure 
4).  Although  the  majority  of  precipitation  at  Fort  Bliss  falls  from  late  summer 
and  early  autumn  thunderstorms,  which  are  often  isolated,  the  El  Paso  precipita¬ 
tion  data  does  provide  a  general  estimate  of  precipitation  at  the  study  site. 
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Figure  2.  Trends  in  spectral  indices  from  MSS  imagery  (1972  through  1984). 
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Figure  4.  Landsat  TM  Tasseled  Cap  brightness  and  annual  precipitation  (1984  through  1995). 
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Although  the  June  1992  scene  was  not  considered  in  evaluating  the  overall  trend 
of  indices  (because  June  varies  from  the  remaining  acquisition  dates),  precipita¬ 
tion  records  indicate  a  relatively  high  amount  of  precipitation  was  recorded  in 
May  1992  (4.22  in.  [10.72  cm]),  which  may  partially  explain  the  significantly 
lower  brightness  indices  observed  for  this  date.  Similarly,  a  relatively  high 
amount  of  precipitation  was  recorded  in  June  (3.05  in.  [7.75  cm])  and  July  (2.66 
in.  [6.76  cm])  of  1986,  again  possibly  explaining  the  observed  decrease  in  bright¬ 
ness  in  the  August  1986  image. 

Conversely,  a  significant  increase  in  brightness  was  observed  for  TCB  and  Red 
Reflectance  in  1987  (Figure  3).  Although  precipitation  for  September  1987  (im¬ 
age  acquisition  date)  was  below  the  mean  precipitation  recorded  in  September 
for  the  time  period  of  this  study  (0.89  in.  [2.26  cm]),  and  precipitation  in  July 
was  also  below  average  (0.64  in.  [1.63  cm]),  total  antecedent  precipitation  leading 
up  to  the  acquisition  date  appeared  to  be  normal  or  even  slightly  above  normal. 
Therefore,  it  does  not  appear  that  this  increase  in  brightness  can  be  totally  ex¬ 
plained  by  existing  precipitation  records;  more  detailed  climatic  data  would  be 
required  to  confirm  this  observation.  It  is  possible  that  the  increase  in  detected 
brightness  in  1987  may  be  the  result  of  increased  training  impact  between  1986 
and  1987,  relative  to  other  temporal  sequences  observed.  However,  cause  and 
effect  requires  more  detailed  field  data,  including  climatic  data,  soil  moisture, 
and  information  on  the  intensity  and  geographic  extent  of  training  exercises 
during  this  time  period.  Note  that  although  spectral  brightness  peaked  in  1987 
according  to  Red  Reflectance  and  TCB,  the  same  increase  was  not  evident  in  Al¬ 
bedo,  which  indicated  a  larger  increase  between  1987  and  1988  (Figure  3). 

In  contrast  to  these  brightness  indices,  a  single  greenness  index  (MSAV12)  was 
also  evaluated.  As  described  earlier,  brightness  indices  are  usually  more  indica¬ 
tive  of  land  condition  than  greenness  indices  in  arid  environments  with  rela¬ 
tively  sparse  vegetation  cover.  Spectral  greenness  and  brightness  indices  are 
usually  inversely  proportional  to  each  other  in  arid  environments.  An  increase 
in  brightness  and  decrease  in  greenness  is  usually  associated  with  degradation, 
while  a  decrease  in  brightness  and  increase  in  greenness  is  usually  associated 
with  an  increase  in  vegetative  cover  and  improvement  in  land  condition.  How¬ 
ever,  in  arid  environments  where  vegetative  cover  does  not  exceed  30  to  35  per¬ 
cent  background  soil  brightness  dominates  the  spectral  response,  and  therefore 
greenness  indices  typically  do  not  correlate  well  with  abundance  of  vegetative 
cover.  Despite  these  recognized  limitations,  MSAVI2  was  evaluated  for  compari¬ 
son. 

Based  on  the  typical  inverse  relationship  between  brightness  and  greenness  in¬ 
dices,  one  would  expect  increases  in  MSAVI2  for  time  periods  where  a  decrease 
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in  brightness  was  observed,  and  vice  versa.  However,  this  only  occurred  between 
1992  to  1993,  and  1993  to  1994,  when  MSAVI  decreased  and  all  three  brightness 
indices  increased,  and  also  in  1987  to  1988,  and  1988  to  1990,  where  MSAVI  in¬ 
creased,  and  some,  but  not  all,  of  the  brightness  indices  decreased  (Figure  3).  In 
June  1992,  when  brightness  values  were  relatively  low,  MSAVI2  values  were 
relatively  high.  However,  in  August  1986,  when  brightness  values  were  rela¬ 
tively  low,  MSAVI2  values  were  also  relatively  low. 

The  minimum  MSAVI2  value  was  recorded  in  1994.  Acquired  at  the  onset  of  se¬ 
nescence,  this  image,  as  well  as  the  November  1992  image,  should  be  interpreted 
with  caution,  or  perhaps  eliminated  from  consideration  when  evaluating  trends. 

Based  on  previous  literature  (Robinove  et  al.  1981;  Musick  1984;  Duncan  et  al. 
1993),  and  the  unexpected  trends  in  MSAVI2  when  compared  with  trends  in 
brightness  indices,  it  was  determined  that  MSAVI2,  as  well  as  other  greenness 
indices  that  were  evaluated  over  the  same  time  period,  were  not  good  indicators 
of  vegetative  cover  and  land  condition.  However,  the  relationship  between 
MSAVI2  and  the  brightness  indices  is  discussed  further  when  land  use  patterns 
are  considered. 

Examination  of  all  three  mean  brightness  indices  from  1972  to  1984  using  MSS 
imagery  reveals  a  similar  pattern,  especially  between  Red  Reflectance  and  TCB 
(Figure  2).  Both  Red  Reflecteince  and  TCB  indices  indicated  that  September 
1972  had  the  highest  reflectance  or  brightness  over  this  time  period,  followed  by 
a  significant  decrease  in  brightness  in  December  1972.  Albedo  was  also  very 
high  in  September  1972,  and  increased  even  higher  by  December  1972,  even 
though  Red  and  TCB  indicated  a  large  decrease  in  brightness  over  this  3-month 
time  period.  Similar  to  Red  and  TCB,  Albedo  decreased  significantly  between 
September  1972  and  October  1974.  Although  there  was  some  continued  slight 
variation  in  the  trend  of  brightness  indices  from  1974  to  1984,  the  September 
1972  image  and,  to  a  lesser  extent,  the  June  1982  image  were  the  only  two  dates 
that  exhibited  relatively  high  brightness  indices.  As  mentioned  earlier,  the  June 
image  should  be  interpreted  with  caution,  as  it  varies  considerably  from  the  tar¬ 
get  acquisition  dates  (August  through  October).  Phenological  differences  be¬ 
tween  the  relatively  dry  period  of  mid-summer  (Jime)  and  the  wet  season 
(August  to  October)  scenes  are  the  most  likely  explanation  for  the  relatively  high 
brightness  indices  observed  for  June  1982. 

Precipitation  records  may  provide  some  explanation  for  the  rapid  decrease  in 
brightness  indices  between  September  1972  and  November  1974.  1974  was  the 
second  wettest  on  record  of  the  23  years  observed  by  this  study  (13.95  in.  [35.43 
cm]).  Also,  September  1974,  which  was  1  month  prior  to  the  image  acquisition 
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date,  was  the  wettest  month  on  record  over  this  same  time  period,  with  6.68 
inches  [16.97  cm]  of  rainfall.  This  amount  of  precipitation  undoubtedly  resulted 
in  a  significant  increase  in  vegetative  cover,  growth,  and  vigor  in  the  following 
month,  as  well  as  an  increase  in  soil  moisture,  both  of  which  would  explain  the 
observed  decrease  in  brightness  during  this  time  period.  The  relatively  high 
amount  of  precipitation  recorded  in  September  1978  may  also  explain  the  rela¬ 
tive  decrease  in  observed  brightness  or  reflectance  in  the  September  1978  image. 

In  contrast  to  these  brightness  indices,  the  greenness  index  (MSAVI2)  increased 
significantly  between  September  1972  and  October  1974.  Again,  this  would  sup¬ 
port  the  hypothesis  that  this  trend  is  most  likely  attributable  to  a  significant  in¬ 
crease  in  precipitation  and  moisture  availability  that  resulted  in  a  substantial 
increase  in  vegetative  cover,  growth,  and  vigor  over  this  time  period. 

Unlike  comparisons  between  trends  in  brightness  and  greenness  indices  derived 
from  TM  imagery,  the  typical  inverse  relationship  between  brightness  and 
greenness  indices  was  evident  for  most  of  the  time  periods  observed  with  MSS 
imagery.  In  addition  to  the  expected  inverse  trends  in  brightness  (decrease)  and 
greenness  (increase)  between  1972  and  1974,  the  trends  reversed  and  were  still 
inverse  to  each  other  between  September  1978  and  June  1982,  when  brightness 
indices  increased  and  greenness  indices  decreases. 


Land  Use  Comparison 

In  addition  to  assessing  relative,  temporal  trends  in  land  condition  using  re¬ 
motely  sensed  spectral  indices  for  the  entire  study  area,  a  second  goal  of  this  re¬ 
search  was  to  compare  trends  for  land  areas  with  contrasting  historical  land 
uses.  As  described  in  Chapter  3,  U.S.  Highway  54  served  as  a  dividing  line  be¬ 
tween  two  areas  of  distinctively  different  historical  land  use  in  the  study  area. 
The  subset  of  the  entire  study  area  that  is  west  of  U.S.  Route  54  (hereafter  called 
West)  is  in  the  Fort  Bliss  Dona  Ana-Orogrande  Complex.  This  subset  includes 
most  of  Maneuver  Areas  5C,  5D,  and  parts  of  4B  and  4C,  all  of  which  have  been 
heavily  impacted  by  tracked  and  wheeled  military  training  and  maneuvering 
from  1972  through  1994.  The  subset  area  to  the  east  of  U.S.  Route  54  (hereafter 
called  East)  is  in  the  McGregor  Guided  Missile  Range.  This  section  has  not  been 
impacted  by  tracked  and  wheeled  military  training  and  maneuvers,  and  there¬ 
fore  served  as  a  control  for  comparison  with  the  heavily  impacted  maneuver  ar¬ 
eas. 

Appendix  B  summarizes  the  same  statistics  for  each  of  the  four  indices  for  the 
East  and  West  subsets  that  were  summarized  for  the  entire  study  area.  In  addi- 
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tion,  mean  index  values  for  East  and  West  for  the  indices  evaluated  for  the  entire 
study  area  were  also  plotted  for  each  image  date  in  Figures  2  and  3. 

Initial  examination  of  all  three  mean  brightness  indices  derived  from  TM  im¬ 
agery  for  the  East  and  West  areas  from  1984  through  1994  revealed  a  trend 
similar  to  that  observed  for  the  full  study  area.  Extremes  in  mean  values  for  the 
East  and  West  land  areas  occurred  during  the  same  years  as  those  observed  for 
the  full  study  area.  Trends  in  mean  brightness  values  for  the  east  and  west  ar¬ 
eas  were  very  similar  to  the  trends  observed  for  the  full  study  area.  This  sug¬ 
gests  that  for  the  years  in  which  extremes  were  observed,  they  were  the  result  of 
climatic  or  phonological  differences  that  occurred  uniformly  across  the  entire 
study  area,  and  were  not  a  result  of  military  impacts.  If  an  increase  in  bright¬ 
ness  for  the  entire  study  area  for  certain  years  was  an  indicator  of  increased 
training  impact,  there  should  have  been  an  observable  difference  in  response  be¬ 
tween  the  impacted  (West)  and  nonimpacted  (East).  This  difference  was  not  ob¬ 
served.  One  observation  of  potential  significance  was  that  the  difference  in 
mean  brightness  between  East  and  West  decreased  substantially  beginning  in 
1992  and  continuing  through  1994.  This  may  indicate  two  possible  trends:  ei¬ 
ther  vegetative  cover  was  increasing  in  the  West  relative  to  the  East,  or  was  de¬ 
creasing  in  the  East  relative  to  the  West.  Since  the  overall  trend  in  brightness 
indices  increased  between  1992  and  1994,  this  would  indicate  that  the  latter  is 
more  likely,  with  vegetative  cover  decreasing  in  both  areas,  but  decreasing  more 
rapidly  in  the  East.  Records  indicated  that  annual  precipitation  was  below  the 
mean  annual  precipitation  for  both  1993  and  1994,  with  1994  being  the  driest 
year  (5.48  in.  [13.92])  of  all  years  observed  in  this  study.  This  may  explain  the 
overall  trend  of  increased  brightness  for  the  entire  study  area,  but  would  not  ex¬ 
plain  the  greater  increase  in  brightness  for  the  East. 

Not  only  was  the  difference  in  brightness  between  East  and  West  relatively  con¬ 
sistent  for  all  years,  with  the  exception  of  1992  through  1994,  but  the  West  area 
had  a  higher  brightness  index  than  the  East  area  for  each  year  observed.  Con¬ 
versely,  the  East  area  had  a  higher  greenness  index  than  the  West  area  for  each 
year  observed.  This  observed  difference  between  East  and  West  corresponds  well 
with  what  was  expected  based  on  the  historical  land  use  of  each  area.  The  West 
has  been  more  heavily  impacted  by  tracked  and  wheeled  vehicles  than  the  East, 
and  therefore,  it  was  expected  that  there  would  be  relatively  less  vegetative 
cover  and  more  exposed  and  disturbed  soil  in  the  West  than  in  the  East.  Conse¬ 
quently,  the  West  area  should  have  a  higher  brightness  index  and  lower  green¬ 
ness  index  than  the  East  area. 

One  final  comparison  of  land  areas  within  the  full  study  area  was  observed  for 
two  areas  that  did  not  have  a  clearly  different  historical  land  use,  but  instead 
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were  visually  identified  on  the  index  images  as  contrasting  areas.  Within  the 
West  area  were  two  visibly  distinct  subareas;  the  first  appeared  more  heavily 
impacted  and  is  referred  to  as  Low  Cover  in  this  study,  the  second  area,  High 
Cover,  was  identified  as  less  impacted. 

Similar  to  comparisons  made  between  East  and  West,  the  same  statistics  for 
these  two  subareas  are  summarized  in  Appendix  B  and  the  annual  mean  index 
values  for  these  areas  are  plotted  in  Figures  2  and  3.  Visually  observed  differ¬ 
ences  between  these  two  areas  were  confirmed  by  assessing  differences  between 
brightness  values  for  each  area  for  each  year  observed.  The  Low  Cover  area  ex¬ 
hibited  the  highest  brightness  index  of  all  areas  observed  in  this  study,  including 
the  full  study  area,  the  East  and  West  study  areas,  and  the  visually  identified 
High  Cover  area  within  the  West  study  area.  Similarly,  the  Low  Cover  area  also 
had  the  lowest  greenness  index  (MSAVI2)  of  all  areas  observed  in  this  study. 
Examination  of  all  three  mean  brightness  indices  for  the  Low  Cover  and  High 
Cover  areas  from  1984  through  1994  using  TM  imagery  (Figure  3)  revealed  a 
trend  similar  to  that  of  the  full  study  area  and  the  East  and  West  areas.  Ex¬ 
tremes  in  mean  values  for  the  East  and  West  land  areas  occurred  during  the 
same  years  as  those  observed  for  the  full  study  area.  Similar  to  earlier  conclu¬ 
sions,  this  would  suggest  that  observed  extremes  were  probably  the  result  of  cli¬ 
matic  or  phonological  differences  that  occurred  xmiformly  across  the  entire  study 
area  rather  than  to  the  varying  degree  of  military  impacts  that  may  have  oc¬ 
curred  in  the  Low  Cover  and  High  Cover  areas. 


Spatial  Patterns 

One  final  objective  of  this  research  was  to  identify  any  spatial  patterns  of  small 
scale,  gross  level  change  that  were  evident  in  the  temporal,  archival  multispec- 
tral  image  data  set.  Several  commonly  used  change  detection  methods  have 
been  applied  to  temporal  image  data  sets  with  varying  results,  as  described  in 
Chapter  2.  The  most  commonly  applied  method,  and  the  one  used  in  this  re¬ 
search,  is  referred  to  as  image  differencing. 

As  the  name  suggests,  images  from  two  different  dates  are  spatially  coregistered 
to  each  other  and  subtracted  using  image  algebra.  The  difference  between  pixel 
values  between  the  two  dates  in  a  temporal  sequence  is  then  assigned  as  the 
pixel  value  in  the  difference  image.  If  the  intent  is  to  monitor  some  specific  as¬ 
pect  of  the  Earth’s  surface,  such  as  vegetative  cover  and/or  exposed  soil,  it  is  ad¬ 
vantageous  to  first  calculate  an  index  from  the  two  original  image  dates  that 
provides  a  good  indicator  of  the  variable  of  interest  on  the  surface.  In  this  case, 
the  index  images  are  subtracted,  rather  than  the  digital  numbers  (DNs)  in  the 
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original  images.  At  a  minimum,  the  DNs  should  always  be  calibrated  to  account 
for  possible  atmospheric  and  radiometric  differences  between  image  acquisition 
dates.  In  many  cases,  a  threshold  is  established  that  represents  “normal” 
change  in  indices  between  two  dates.  Areas  that  increase  or  decrease  beyond 
this  threshold  are  identified  as  improved  or  degraded.  Establishing  a  threshold 
allows  identification  of  areas  that  have  increased  or  decreased  beyond  the  mean 
change  between  two  image  dates.  This  serves  to  eliminate  changes  that  occurred 
uniformly  across  an  image,  such  as  climatic  or  phonological  change,  and  isolates 
areas  where  there  has  been  an  actual  change  in  land  condition.  Selection  of 
similar  image  acquisition  dates  for  different  years  helps  to  eliminate  or  minimize 
phonological  differences  as  well. 

Using  image  index  differencing,  patterns  of  change  were  identified  for  three  dif¬ 
ferent  temporal  sequences  using  TM  imagery.  Brightness  indices,  and  not 
Greenness  indices,  were  selected  for  image  differencing  based  on  previous  re¬ 
search  and  literature  (Robinove  et  al.  1981;  Musick  1984;  Duncan  et  al.  1993), 
which  indicates  that  brightness  indices  are  best  for  monitoring  vegetative  cover 
in  arid  environments.  Differences  between  several  brightness  indices  were 
mapped  and  evaluated,  and  it  was  determined  that  similar  patterns  were  de¬ 
tected  using  different  brightness  indices.  Red  Reflectance  and  TCB,  were  se¬ 
lected  for  detailed  analysis  and  presentation.  Changes  in  these  two  indices  were 
mapped  between  1987  to  1990,  1990  to  1993,  and  overall  from  1987  to  1993.  For 
each  time  sequence,  the  mean  difference  between  the  two  image  dates  was  de¬ 
termined,  and  areas  that  increased  or  decreased  beyond  2  standard  deviations 
from  the  mean  change  were  mapped  within  the  study  area.  These  images  are 
presented  in  Figures  5,  6,  and  7. 

As  mentioned  earlier,  an  increase  in  a  brightness  index  is  usually  associated 
with  degradation  or  loss  of  vegetative  cover,  while  a  decrease  in  a  brightness  in¬ 
dex  is  usually  associated  with  an  increase  in  vegetative  cover.  Therefore,  those 
areas  that  decreased  by  more  than  2  standard  deviations  from  the  mean  de¬ 
crease  between  the  two  dates,  or  those  areas  that  may  have  increased,  but  by 
less  than  2  standard  deviations  from  the  mean  increase  were  highlighted  in 
green  on  the  map  to  signify  a  potential  increase  in  vegetative  cover  and  im¬ 
provement  in  land  condition.  Conversely,  those  areas  that  increased  more  than  2 
standard  deviations  from  the  mean  increase,  or  those  areas  that  may  have  de¬ 
creased,  but  by  less  than  2  standard  deviations  from  the  mean  decrease  were 
highlighted  in  red  on  the  map  to  signify  a  potential  decrease  in  vegetative  cover 
and  degradation.  However,  these  areas  of  delineated  change  must  be  interpreted 
with  caution,  as  there  was  no  field  data  to  verify  and  validate  the  observed 
changes.  Therefore,  cause  and  effect  was  difficult,  if  not  impossible,  to  establish. 
Several  other  factors  that  are  not  related  to  actual  change  in  land  condition  may 
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to  these  areas  of  observed  change.  Simple  spatial  registration  or  atmospheric 
calibration  errors  between  image  dates  could  result  in  the  misclassification  of 
change.  Spatial  registration  errors  are  most  noticeable  along  linear  features 
such  as  roads  or  other  features  with  a  high  contrast  edge  between  two  land  cover 
types.  In  addition,  even  though  images  of  similar  anniversary  dates  were  ana¬ 
lyzed,  there  may  still  have  been  phenological  differences  between  the  dates.  And, 
there  may  still  have  been  differences  in  climatic  conditions  between  the  two  sea¬ 
sons,  including  differences  in  precipitation  in  the  preceding  months,  or  even  a 
very  recent  precipitation  event,  which  would  have  resulted  in  increased  soil 
moisture  for  one  of  the  dates. 

Between  1987  and  1990,  there  was  a  decrease  in  mean  brightness  for  both  TCB  (- 
4.83)  and  Red  Reflectance  (-23.41)  for  the  entire  study  area.  Figure  5  shows  a 
map  of  those  areas  of  extreme  change  (+/-  2  s.d.  from  the  mean  change)  for  TCB 
and  Red  Reflectance.  Green  highlighting  identifies  those  areas  that  became 
much  less  bright  or  less  reflective  over  this  time  period;  conversely,  red  highlight¬ 
ing  identifies  those  areas  that  became  much  brighter  or  more  reflective  over  this 
same  time  period.  Patterns  of  extreme  change  were  very  similar  for  both  bright¬ 
ness  indices,  which  was  expected.  Several  areas  of  observed  extreme  change 
were  visited  in  the  field.  Although  field  verification  was  conducted  several  years 
after  the  time  of  the  observed  change,  field  observations  confirmed  that,  in  most 
cases,  areas  of  extreme  change  seemed  to  occur  in  depressions,  or  isolated  areas 
of  high  productivity  and  vegetative  cover,  or  in  areas  of  extremely  low  vegetative 
cover  or  completely  exposed  soils. 

Between  1987  and  1990,  those  areas  that  decreased  much  less  than  the  mean 
decrease  in  brightness  for  the  study  area  (highlighted  in  green)  appeared  to  oc¬ 
cur  in  isolated  areas  of  extremely  low  vegetative  cover,  or  bare  ground.  This 
would  generally  indicate  an  increase  in  vegetative  cover  or  decrease  in  exposed 
soils  for  these  areas.  Precipitation  records  indicate  that  precipitation  before  the 
1990  date  was  significantly  higher  than  precipitation  before  the  1987  date. 
Therefore,  these  areas  would  have  had  more  moisture  available,  and  as  a  result, 
one  would  expect  an  increase  in  vegetative  cover,  growth,  and  vigor  in  these  ar¬ 
eas  in  1990  relative  to  1987.  It  is  likely  that  these  areas  experienced  a  combina¬ 
tion  of  green-up  of  perennials  and  annuals,  an  increase  in  soil  moisture,  and  pos¬ 
sibly  an  increase  in  shadows,  which  may  have  caused  these  areas  to  appear 
significantly  darker  in  1990.  Areas  where  these  extreme  changes  were  most  evi¬ 
dent  were  on  the  McGregor  Range  and  directly  west  of  Highway  54  in  Maneuver 
Area  5D,  as  well  as  a  few  isolated  areas  of  extremely  low  cover  in  Maneuver  Area 
5C.  These  areas  also  occurred  along  many  of  the  roads  and  trails.  There  may 
have  been  an  increase  in  soil  moisture  in  these  areas,  or  the  extreme  change  may 
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During  this  same  time  period,  the  areas  that  increased  in  brightness  or  de¬ 
creased  much  less  than  the  mean  decrease  in  brightness  for  the  study  area  (high¬ 
lighted  in  red)  appeared  to  occur  in  isolated  areas  of  high  productivity  or  vegeta¬ 
tive  cover.  An  increase  in  brightness,  or  a  decrease  much  less  than  the  mean 
decrease  would  generally  indicate  a  decrease  in  vegetative  cover  or  increase  in 
exposed  soils  for  these  areas,  but  because  these  areas  of  extreme  change  occurred 
in  areas  of  high  cover,  it  is  more  likely  that  these  extreme  changes  were  due  to 
differences  in  climatic  conditions  between  the  two  dates.  However,  precipitation 
records  indicate  that  precipitation  before  the  1990  date  was  significantly  higher 
than  precipitation  before  the  1987  date.  Therefore,  these  areas  of  high  productiv¬ 
ity  should  have  had  more  moisture  available,  especially  in  depression  areas,  and 
one  would  expect  an  increase,  and  not  a  decrease,  in  vegetative  cover  and  vigor  in 
these  areas.  It  is  probably  more  likely  that  a  single  precipitation  event  just  be¬ 
fore  the  1987  image  acquisition  date  caused  a  green-up  of  perennials  and  annu¬ 
als  and  an  increase  in  soil  moisture,  which  may  have  caused  these  areas  to  ap¬ 
pear  significantly  brighter  in  1990.  Areas  where  these  extreme  changes  were 
most  evident  were  the  extreme  western  edge  of  the  study  area,  extending  south 
from  Old  Coe  Lake  in  Maneuver  Area  4B,  and  to  a  lesser  extent,  in  Maneuver 
Area  4C.  A  second  area  where  this  type  of  change  was  evident  was  in  the  south¬ 
east  comer  of  the  study  area,  on  the  McGregor  Range,  directly  east  of  Alvarado 
Tank.  Other  similar  changes  occurred  in  isolated  areas  of  high  cover  throughout 
the  study  area. 

Between  1990  and  1993,  there  was  a  very  slight  increase  in  mean  brightness  for 
both  TCB  (0.10)  and  Red  Reflectance  (10.84)  for  the  entire  study  area.  Figure  6 
is  a  map  of  those  areas  of  extreme  change  (+/-  2  s.d.  from  mean  change)  for  TCB 
and  Red  Reflectance.  Similar  to  the  1987  to  1990  evaluation,  areas  of  extreme 
change  occurred  in  primarily  the  same  areas.  However,  the  same  areas  that  be¬ 
came  less  bright  or  less  reflective  between  1987  and  1990  seemed  to  increase  in 
brightness  and  reflectivity  between  1990  and  1993.  Likewise,  those  areas  that 
increased  in  brightness  or  reflectivity  between  1987  and  1990  seemed  to  decrease 
in  brightness  and  reflectivity  between  1990  and  1993.  Besides  this  reversal  in 
direction  of  change,  the  only  other  notable  difference  between  the  two  different 
time  sequences  appeared  to  be  an  increase  in  total  area  of  extreme  change  on  the 
McGregor  Range  from  1990  to  1993  relative  to  1987  to  1990.  Unlike  1987  to 
1990,  where  there  appeared  to  be  a  relatively  equal  amount  of  extreme  increases 
and  decreases,  between  1990  and  1993,  almost  all  of  the  extreme  changes  were 
increases  in  brightness.  These  extreme  increases  in  brightness  indicate  a  rela¬ 
tive  decrease  in  vegetative  cover  or  vigor,  which  may  have  been  the  result  of  less 
antecedent  precipitation  before  the  1993  acquisition  date  relative  to  the  1990  ac¬ 
quisition  date. 


Legend : 


Increase  less  than  2  standard  deviations  from  mean  change  or  decrease 
2  in  Landsat  TM  Band  3  (Red)  Reflectance  between  1990  and  1993. 

Increase  in  Landsat  TM  Band  3  (Red)  Reflectance  greater  than 
2  standard  deviations  from  mean  change  between  1990  and  1993. 


Figure  6.  Spatial  patterns  of  extreme  change  in  Red  Reflectance  and  TCB  (1990  through  1993). 

Between  1987  and  1993,  there  was  a  decrease  in  mean  brightness  for  both  TCB  (- 
4.72)  and  Red  Reflectance  (-12.57)  for  the  entire  study  area.  Figure  7  shows  a 
map  of  those  areas  of  extreme  change  {+!-  2  s.d.  from  mean  change)  for  TCB  and 
Red  Reflectance.  Spatial  patterns  of  extreme  change  for  this  time  period  were 
very  similar  to  those  observed  between  1990  and  1993.  Again,  most  of  the  areas 
of  extreme  change  were  areas  that  either  increased  in  brightness,  or  decreased, 
but  much  less  than  the  mean  decrease  for  the  entire  study  area.  Most  of  these 
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increases  in  brightness  occurred  in  north-south  oriented  linear  depressions  that 
extend  south  from  Old  Coe  Lake  in  Maneuver  Area  4B,  and  in  the  southeastern 
comer  of  the  study  area,  on  McGregor  Range,  directly  east  of  Alvarado  Tank.  In 
each  of  the  time  periods  observed,  there  appeared  to  be  an  area  of  increasing 
brightness,  or  decreasing  cover,  in  the  southeastern  comer  of  the  study  area,  on 
McGregor  Range.  Between  1987  and  1993,  the  total  area  of  increased  brightness 
was  the  greatest  for  all  time  periods  observed.  Since  this  area  of  McGregor 
Range  is  not  impacted  by  tracked  and  wheeled  vehicles,  it  is  more  likely  that  the 
observed  patterns  of  increasing  brightness  are  due  to  natural  changes  rather 
than  military  impacts. 

The  spatial  patterns  of  extreme  change  between  different  dates  were  identified 
by  delineating  those  areas  that  changed  beyond  2  standard  deviations  from 
mean  change.  By  evaluating  differences  from  the  mean  change,  rather  than  the 
simple  difference  between  two  dates,  those  climatic  effects  that  would  have  af¬ 
fected  brightness  values  uniformly  across  the  entire  study  area  should  have  been 
minimized.  However,  the  extreme  changes  that  were  observed  were  still  most 
likely  attributable  to  climatic  factors,  such  as  differences  in  soil  moisture,  or  dif¬ 
ferences  in  vegetative  growth,  vigor,  and  greenness.  The  differences  were  simply 
more  extreme  in  isolated  areas  of  either  high  or  low  vegetative  cover. 


Decrease  in  Tasseled-Cap  Brightness  Index  greater  than 
2  standard  deviations  from  mean  change  between  1987  and  1993. 

Increase  in  Tassled-Cap  Brightness  Index  greater  than 
2  standard  deviations  from  mean  change  between  1987  and  1993. 


Figure  7.  Spatial  patterns  of  extreme  change  in  Red  Reflectance  and  TCB  (1987  through  1993). 


ERDC/CERLTR-01-36 


45 


6  Summary  and  Recommendations 


Summary 

The  primary  goal  of  this  research  was  to  assess  relative,  temporal  trends  in  land 
condition  using  remotely  sensed  spectral  indices  for  a  specific  study  site  at  Fort 
Bliss,  Texas.  A  secondary  goal  of  this  research  was  to  observe  spatial  patterns  of 
detected  changes  in  land  condition,  and  to  compare  trends  and  patterns  for  areas 
that  have  been  impacted  by  military  training  with  those  that  have  not  been  im¬ 
pacted  by  the  training  mission. 

Spectral  brightness  and  greenness  indices  were  calculated  from  temporal  im¬ 
agery  and  were  used  as  a  surrogate  measure  of  the  relative  abundance  of  vegeta¬ 
tive  cover  on  the  Earth’s  surface.  Observations  of  trends  in  brightness  and 
greenness  indices  over  a  23-year  period  from  1972  through  1994  revealed  consid¬ 
erable  variations  in  brightness  or  reflectance  of  the  study  area.  However,  when 
analyzed  in  conjunction  with  precipitation  data  for  the  same  time  period,  most 
variation  in  spectral  indices  could  be  related  to  natural  variations  in  the  amount 
of  monthly  and  annual  precipitation  for  the  region.  In  most  cases,  large  in¬ 
creases  or  decreases  in  spectral  indices  were  probably  the  result  of  an  increase  in 
growth  and  vigor  of  both  perennials  and  annuals  during  seasons  of  above  normal 
precipitation,  and  conversely,  a  decrease  in  growth  and  vigor  during  drought  pe¬ 
riods.  In  addition,  although  extremes  in  brightness  and  greenness  indices  were 
observed  for  some  individual  image  dates,  spectral  indices  appeared  to  return  to 
close  to  mean  values  after  a  relatively  short  period  of  time.  This  would  seem  to 
indicate  that  land  condition  tended  to  degrade  or  improve  in  response  to  short¬ 
term  climate  variations,  or  possibly  anthropogenic  disturbances  such  as  military 
training,  but  there  was  no  indication  of  permanent,  long-term  changes  in  land 
condition.  In  summary,  the  study  area  appeared  to  be  resilient,  with  no  clearly 
observable  trends  toward  either  long-term  degradation  or  improvement. 

Comparison  of  trends  in  spectral  indices  between  impacted  and  unimpacted  ar¬ 
eas  of  the  study  site  revealed  a  noticeable  difference  between  the  heavily  im¬ 
pacted  maneuver  areas  on  the  west  side  of  Highway  54  and  the  McGregor 
Guided  Missile  Range  on  the  east  side  of  the  highway.  Brightness  indices  for  the 
maneuver  areas  were  consistently  higher  than  brightness  indices  for  McGregor 
Range,  indicating  that  the  maneuver  areas  most  likely  have  less  total  vegetative 
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cover.  However,  the  observed  trends  in  brightness  and  greenness  indices  for  both 
the  East  and  West  study  areas  were  very  similar.  This  indicated  that  the  cause 
of  these  observed  variations  in  trend  was  more  likely  attributable  to  climatic  or 
phonological  changes  that  occurred  uniformly  across  the  study  area  as  opposed 
to  military  or  other  anthropogenic  impact.  If  an  increase  in  brightness  provided 
some  indication  of  a  decrease  in  vegetative  cover  or  increase  in  soil  disturbance 
due  to  military  impacts,  there  should  have  been  an  observable  difference  in  re¬ 
sponse  between  the  maneuver  areas  (West)  and  McGregor  Range  (East).  This 
difference  was  not  observed.  The  only  exception  to  this  observed  trend  was  be¬ 
tween  1992  and  1994,  where  differences  between  indices  for  the  East  and  West 
decreased  substantially. 

Observations  during  the  time  period  of  this  study  do  not  provide  enough  data  to 
conclude  that  observed  differences  between  East  and  West  can  be  attributed  to 
military  impact.  The  contrasting  historical  land  uses  between  the  maneuver  ar¬ 
eas  and  McGregor  Range  provide  the  most  reasonable  explanation  as  to  why 
there  is  a  considerable  difference  in  spectral  indices  for  the  two  areas.  However, 
during  the  time  period  of  this  study,  there  has  not  been  a  significant  increase  or 
decrease  in  the  difference  between  these  two  areas.  The  differences  are  probably 
the  result  of  cumulative  training  impacts  over  a  much  longer  time  period.  Com¬ 
parison  of  trends  in  indices  between  an  area  of  high  and  low  vegetative  cover 
within  the  maneuver  areas  revealed  similar  results.  Again,  there  was  a  consis¬ 
tent  difference  in  spectral  indices  for  the  two  areas,  but  no  significant  variation 
in  the  difference  between  the  two  areas  that  would  indicate  that  a  change  had 
actually  occurred  during  the  time  period  of  this  study. 

Spatial  patterns  of  change  were  variable  across  different  temporal  sequences,  as 
expected.  Different  change  detection  techniques,  when  applied  to  the  same  tem¬ 
poral  sequence,  revealed  similar  patterns  of  change,  which  served  to  verify  that 
detected  changes  in  the  imagery  represented  actual  changes  in  land  condition 
and  cover  on  the  ground.  In  most  cases,  areas  of  extreme  change  seemed  to  oc¬ 
cur  in  depressions,  or  isolated  areas  of  high  productivity  and  vegetative  cover,  or 
in  areas  of  extremely  low  vegetative  cover  or  completely  exposed  soils.  There 
was  no  indication  that  these  areas  of  extreme  change  occurred  in  areas  that  were 
heavily  impacted  by  tracked  and  wheeled  vehicle  maneuvers.  Again,  the  ex¬ 
treme  changes  in  spectral  indices  that  were  observed  for  these  areas  were  most 
likely  attributable  to  short-term  changes  in  vegetative  growth  and  vigor  in  re¬ 
sponse  to  short-term  variations  in  the  amount  of  rainfall  and  corresponding  soil 
moisture. 

Quantification  of  the  abundance  of  vegetative  cover  with  remotely  sensed  im¬ 
agery  is  difficult  and  requires  an  empirical  relationship  between  indices  and 
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ground  measurements  for  the  purpose  of  calibration.  Field  measurements  of 
cover  were  not  available  for  calibration;  therefore,  only  relative  trends,  rather 
than  absolute  changes  in  cover  amounts,  were  assessed  in  this  study.  Even  if 
ground  observations  of  vegetative  cover  had  been  available,  the  archival  imagery 
from  which  the  spectral  indices  were  derived  was  collected  at  a  very  coarse  spa¬ 
tial  resolution  (80  m  and  30  m).  In  arid  environments  such  as  Fort  Bliss,  where 
vegetative  cover  is  sparse,  a  single  pixel  records  the  reflectance  of  a  heterogene¬ 
ous  mixture  of  vegetative  cover,  including  both  live  and  standing  dead  material, 
litter,  shadows,  and  significant  areas  of  bare  ground  between  plants.  Therefore, 
because  of  this  mixed  spectral  response,  it  would  have  been  very  difficvdt  to  cali¬ 
brate  spectral  indices  with  ground  measiirements  for  the  purpose  of  calibration 
and  actual  quantification  of  vegetative  abimdance. 

Assessment  of  trends  in  vegetation  and  brightness  indices  calculated  from  tem¬ 
poral  image  data  sets  still  provides  a  cost-effective  method  for  monitoring  rela¬ 
tive  trends  in  land  condition.  Documentation  of  historical  trends  in  land  condi¬ 
tion  may  be  useful  at  an  upper  management  level  for  evaluating  the 
effectiveness  of  historical  resource  management  efforts,  including  Land  Rehabili¬ 
tation  and  Maintenance  over  long  periods  of  time  and  across  large  geographic 
areas.  Documentation  may  also  be  useful  for  identification  of  areas  of  extreme 
degradation,  or  “problem  areas,”  where  there  is  a  continual  trend  toward  degra¬ 
dation  and  no  evidence  of  recovery.  The  ability  to  monitor  vegetative  cover  dy¬ 
namics  of  training  lands  may  provide  more  information  about  land  condition  and 
carrying  capacity  of  these  lands  than  actual  measurements  of  vegetative  cover, 
structure,  or  pattern  at  any  one  time.  Therefore,  although  the  trends  observed 
in  this  study  provided  only  a  relative,  rather  than  absolute,  measure  of  changing 
vegetative  cover,  they  are  still  of  great  utility  to  installation  land  managers. 


Recommendations 

Archival  remotely  sensed  imagery,  and  specifically,  soil  brightness  and  vegeta¬ 
tion  indices  calculated  from  remotely  sensed  reflectance,  provide  installation 
land  managers  with  the  necessary  tools  to  assess  relative  changes  in  general 
land  condition  at  small  scales,  or  across  large  geographic  areas.  Characteriza¬ 
tion  and  monitoring  of  absolute  changes  in  vegetative  cover  requires  detailed 
ecological  field  surveys  that  allow  for  empirical  relationships  between  remotely 
sensed  reflectance  and  ground-based  cover  measures.  Alternatively,  the  use  of 
spectral  mixture  models,  which  estimate  the  amount  of  cover  based  on  a  propor¬ 
tional  amount  of  the  total  reflectance  that  originates  from  vegetative  cover,  could 
provide  improved  estimates  of  absolute  change  in  vegetative  cover.  Higher  spa¬ 
tial  resolution  imagery  is  also  optimal  for  characterization  of  cover  in  arid  envi- 
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ronments  where  vegetative  cover  is  sparse  and  variable  over  very  short  dis¬ 
tances. 

Regardless  of  whether  the  intent  is  to  monitor  relative  or  absolute  changes  in 
vegetative  cover  on  military  lands,  establishing  cause  and  effect  relationships 
requires  not  only  detailed  ecological  field  surveys,  climate  data,  and  higher  spa¬ 
tial  and  spectral  resolution  imagery,  but  also  detailed  information  on  military 
training  activities.  Improved  information  on  the  spatial  and  temporal  distribu¬ 
tion  and  intensity  of  training  and  testing  activities  is  required  to  establish  rela¬ 
tionships  between  military  impact  and  observed  patterns  of  change  in  land  con¬ 
dition  using  combined  remote  sensing  and  field  survey  technologies.  Collection 
of  ecological  data,  as  well  as  training  data,  is  both  costly  and  time  consuming, 
especially  for  large  installations  such  as  Fort  Bliss.  Therefore,  methods  are  re¬ 
quired  to  scale  up  absolute  estimates  of  vegetative  cover  derived  from  spectral 
mixture  models  and  cover-radiance  models  to  smaller  scales  or  larger  geographic 
regions  for  monitoring  purposes.  In  the  absence  of  reliable  procedures  for  scaling 
between  large  scale  observations  at  field  sample  sites  and  small  scale  observa¬ 
tions  of  entire  training  areas  or  installations,  temporal  analysis  of  coarse  resolu¬ 
tion  satellite  imagery  still  provides  a  cost-effective  method  for  relative  charac¬ 
terization  and  monitoring  of  installation  training  lands. 
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Appendix  A:  Calibration  Values 

Digital  Numbers  extracted  from  a  pond  in  Sunlake  Park,  El  Paso,  Texas. 


MSS  Images 


Acq.  Date: 

Band  1 

Band  2 

Band  3 

Band  4 

09-06-72 

13.906 

10.171 

8.588 

13.221 

12-23-72 

17.546 

15.868 

15.160 

13.883 

10-20-74 

15.011 

15.962 

13.527 

9.605 

11-25-74 

16.821 

14.787 

14.906 

13.120 

09-30-76 

12.735 

13.457 

15.688 

15.025 

09-29-78 

12.015 

12.604 

13.008 

12.983 

12-26-80 

11.983 

11.009 

10.862 

8.310 

06-28-82 

11.451 

*7.983 

*7.129 

*6.658 

12-05-83 

*10.559 

9.587 

9.165 

9.033 

10-28-84 

13.491 

11.484 

10.581 

9.474 

TM  Images 


Acq.  Date: 

Band  1 

Band  2 

Band  3 

Band  4 

Band  5 

Band  7 

10-28-84 

10.97 

8.627 

7.174 

4.272 

0.918 

*0.086 

08-15-86 

*5.310 

*3.358 

*2.401 

*1.443 

*0.734 

0.197 

09-19-87 

9.900 

7.621 

6.256 

4.342 

1.938 

1.149 

08-24-88 

6.760 

4.305 

3.108 

2.024 

1.123 

0.350 

10-13-90 

6.502 

6.051 

4.338 

2.099 

0.866 

0.428 

06-04-92 

7.597 

6.757 

4.965 

4.873 

2.790 

1.917 

10-18-92 

6.630 

6.170 

4.424 

2.141 

0.884 

0.436 

09-19-93 

5.959 

4.390 

3.567 

2.774 

1.544 

0.958 

11-09-94 

6.437 

4.544 

3.030 

1.733 

0.967 

0.582 

*  an  asterisk  indicates  lowest  mean  value  for  each  band  to  subtract  from  other  images. 
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